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. 5^ i As a statistical measure to quantify the topological structure of the large-scale struc- 

ture in the universe, the genus number is calculated for a number of non-Gaussian dis- 
tributions in which the density field is characterized by a nontrivial function of some 
Gaussian-distributed random numbers. As a specific example, the formulae for the log- 
normal and the chi-square distributions are derived and compared with the results of 
A'^-body simulations together with the previously known formulae for the Gaussian dis- 
tribution and second-order perturbation theory. It is shown that the lognormal formula 
fits most of the simulation data the best. 
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1 Introduction 



One of the most important purposes of observational cosmology is to probe the power 
spectrum and the statistical distribution of primordial fluctuations out of the redshift 
survey of large-scale structure, in order to single out the correct model of the evolution 
of our universe. A number of measures have been used to characterize statistical prop- 
erties of the large-scale structures. The most commonly used quantity is the two-point 
correlation function (Totsuji & Kihara 1969), which is nothing but the Fourier trans- 
form of the power spectrum and does not contain much information on its statistical 
distribution. On the other hand, the counts-in-cells statistics, including the void proba- 
bility (White 1979), directly measure the statistical distribution of galaxies and various 
theoretical models based on physical or mathematical arguments have been proposed to 
fit the observational data (Fry 1986). Although this measure contains mathematically 
full information of the statistics in principle, it is difficult to relate the count analysis 
with the visual image or connectivity of galaxy clustering such as filamentary networks, 
sheet-like or bubble- like structures, etc.. 

As a statistical measure to characterize such a topological structure of galaxy distri- 
bution, the genus number has been widely used in the analysis of recent redshift surveys 
(Gott, Melott & Dickinson 1986; Gott, Weinberg & Melott 1987; Weinberg, Gott & 
Melott 1987; Melott, Weinberg & Gott 1988; Gott et al. 1989; Park & Gott 1991; Park, 
Gott & da Costa 1992; Weinberg & Cole 1992; Moore et al. 1992; Vogeley et al. 1994; 
Rhoads, Gott &; Postman 1994). Theoretically, however, the value of the genus had been 
calculated only for the random Gaussian field (Adler 1981; Doroshkevich 1970; Bardeen 
et al. 1986; Hamilton et al. 1986) for a long time except for the restricted cases such 
as Rayleigh-Levy random-walk fractal (Hamilton 1988) and union of overlapping balls 
(Okun 1990). Hence what one could have done at best was to estimate the genus num- 
ber on large scales, which are still in the linear regime, to test the validity of random 
Gaussian initial conditions. 

The situation was somewhat improved recently because lowest-order correction to 
the Gaussian genus number was analytically obtained by one of us using the multi- 
dimensional Edgcworth expansion around the Gaussian distribution (Matsubara 1994). 
Detailed comparison has also been done with the results of A^-body simulations and it 
has been shown that the new formula fits the numerical data well in the semi-linear 
regime but not in the nonlinear regime (Matsubara & Suto 1995). This is in accord 
with the fact that the one-point probability distribution function (PDF) based on the 
Edgeworth series ceases to fit the counts-in-cells once the root-mean-square (rms) value 
of the density contrast becomes as large as ~ 1/4 (Juszkiewicz et al. 1995; Ueda & 
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Yokoyama 1995). Thus it is also desirable to find an analytic expression of the genus 
number for realistic non-Gaussian distributions just as various non-Gaussian models have 
been proposed for counts- in-cells statistics (Fry 1986). 

In the present paper we present analytic formulae of genus number for some non- 
Gaussian distributions in which the density contrast is given by a nontrivial function 
of either a single Gaussian-distributed random number or its combinations. The former 
includes the lognormal distribution which has been used as a model of galaxy distribution 
ever since Hubble (1934). The lognormal model is found to fit the three-dimensional PDF 
of observed galaxies (Hamilton 1985; Bouchet et al. 1993; Kofman et al. 1994), and of 
cell-counts in CDM-type A^-body simulations (Kofman et al. 1994; Ueda & Yokoyama 
1995). Coles & Jones (1991) argued for the lognormal mapping of the linear density 
field to describe its nonlinear evolution in the universe. The latter includes chi-square 
distribution which is closely related with the negative binomial distribution, another 
widely used distribution with a hierarchical property of higher-order cumulants (Fry 
1986; Carruthers 1991; Gaztahaga & Yokoyama 1993; Bouchet et al. 1993). We then 
compare the new formulae with A^-body simulation data with various initial spectra as 
well as the previously known formulae for the Gaussian distribution and second-order 
perturbation theory. 

The rest of the paper is organized as follows. In §^ we review derivation of the 
genus number in Gaussian distribution and its lowest-order correction arising from the 
three-point correlation function. In §^ the genus number is given in the case the density 
contrast is given by a nontrivial but monotonic function of a Gaussian-distributed quan- 
tity. Then in ^ it is extended to the case the density fluctuation is a function of several 
independent Gaussian variables. In §|5| various analytic formulae of the genus and the 
PDF are compared with numerical results obtained from A^-body simulations. Finally 
§^ is devoted to discussion and conclusions. 

2 Genus Curve in Gaussian Distribution and its Nonlinear 
Correction through the Edgeworth Series 

The genus curve G^v) is defined by —1/2 times the total Euler number per unit volume 
of isodensity contours of a continuous density field. It corresponds to 

(holes) — # (isolated regions)] /volume. 

Here the density threshold of the contour is specified hy u = bj \J (S^) with b being the 
density contrast which is smoothed appropriately. Mathematically, it is given by the 
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following expectation value (Doroshkevich 1970; Adler 1981; Bardeen et al. 1986). 

G(^) = -\ {6d (5(x) - i^a) <5D(r?i)'^D(%)|r?3|(CiiC22 - Cu)) , (1) 

where r]i = di5{x.) = 5i(x), Qj = didj6{x) = 5ij(x) and 5d(') is Dirac's delta function. 
Thus it can in principle be calculated once a seven-point probability distribution function 
of 5{x) is known. 

In the case (5(x) obeys Gaussian distribution, we find 6{x), rjiix), and Cjki^) are also 
Gaussian-distributed with a vanishing mean and two-body correlations given by 



(^^(x)) = a', (J(x)^.(x))=0, (5(x)0,(x)) = -^<):- 
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(r?i(x)77,(x)) = ^6ij, (r?i(x)Ofc(x)) = 0, (2) 

(Oj(x)CA:z(x)) = -^{Sijhl + 5ik5jl + Sil6jk), 

where a^, af and a"^ are defined, respectively, by 

a^^{6'), al^{iV6f), and al^UVHf). (3) 

The final result is 

'^M = I^F (1 - ''^) - GkcM. (4) 

If the initial fluctuation is a Gaussian random field, linear theory predicts that the 
genus is described by this Gaussian formula (§) and the PDF by 

e-^V2 



Next we consider correction due to the presence of higher-order correlation func- 
tions, which arises as a result of nonlinear gravitational evolution, in terms of the multi- 
dimensional Edgeworth expansion following Matsubara (1994) and Matsubara & Suto 
(1995). The final result is 



— J^ / 2 \ 2 

e 2 / ar \ 



(27r)2 \ 3cj2 



(6) 

In the above expression, Hn{y) = (— )"e'^ /^{d/dh')'^e~'^ tj^e n-th order Hermite 
polynomial, and S, T, and ?7 are defined as 

1 



1 



T = -^(.2V2.), (7) 
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U = 777w [7(3 + 2i^)^440 + 2IL351 - 5(3 + 4K)L442 - 2IL353 - 6(1 - K)Lu4] 



respectively, which we cah generahzed skewness. They come from the three-point corre- 
lation function that can be evaluated also by second-order perturbation theory (Matsub- 
ara 1994). One should note that generalized skewness should be evaluated for smoothed 
fluctuations that are evolved nonlinearly. Only after taking into account such smooth- 
ing effect, one can compare the second-order perturbation theory with observations. If 
one uses the Gaussian window with the smoothing length R, and assumes the Gaussian 
initial fluctuations, they are explicitly computed as 

•5 = 4^ [(2 + K)L220 + 3Li3i + (1 - K)L222] , 

T = g^[5(5 + 2i^)L240 + 3(9 + K)L33i + 15Li5i 

+ 10(2 -/^)L242 + 3(1 -i^)L333], (8) 

1 

14"0^4 

Here Lai3n{R) stands for the following integral 

L^pn{R)^'^z^ dxj^ dyJ_^df,e-'''^-'+y'+^^-'y^x"y^PMP{x)P{y) (9) 

(10) 

where a and ai are defined by (^) in which 6 is the Gaussian smoothed density fluctuation 
of linear theory over the scale R, Pn is the n-th order Legendre polynomial, and ly is a 
modified Bessel function. The above results (^) to (^) hold for arbitrary values of density 
parameter Q. and cosmological constant A. The latter effect manifests only through the 
function K = K{^}, A) which very weakly depends on and A (Bouchet et al. 1992; 
Bernardeau 1994) where A = A/(3-ff^), H is the Hubble parameter. Its explicit form for 
K has been derived by Matsubara (1995) as 

K{n, A) = ^ - ^ - ^ 3 ^1 ^^^-3/2^ J^' dxX-'/\ (11) 

where 

x{x) = n/x + + X. (12) 

In the two specific models we adopt below, we find if (1,0) = 3/7 = 0.4286 and 
if(0.2,0.8) = 0.4335. 

For the power-law fluctuation spectra P{k) oc k"^, S, T and U can be written down 
explicitly in terms of the hyper geometric function as 



2'4/ ' ' V 2 ' 2 '2'4 



5 



3F 



n 



_3 n + 5 3 1 



+- 



n 



- {n + ?,-K)F 
2){l-K) 



n + 3 n + 5 5 1 
" '2'4 



U 



n+5 n+5 5 1 
2 ' 2 '2'4 



15 
n + 4 



4i^ 



-F 



2 ' 2 
n + 3 n + 5 7 1 
2 ' 2 '2'4 
n + 5 n + 5 7 1 
2 ' 2 '2' 4 



(13) 



The expressions for S in equations and (|T^) are derived by Lokas et al. (1994) which 
are equivalent to the other form independently derived by Matsubara (1994). Similarly 
we transform the expressions for T and U presented in Matsubara (1994; eqs. [16] and 
[18]) using the function La^niR), which are given in equations (^) and (p^). 

This result of equation (|6|) is the analog of the second-order Edgeworth series of PDF 
(Juszkiewicz et al. 1995; Bernardeau & Kofman 1995): 



-1/2/2 



27r 







(14) 



3 The Case Density Field is a Monotonic Function of a 
Random Gaussian Variable 

Next we turn to evaluation of G(z^) for a non-Gaussian distribution whose statistical 
properties are characterized by a monotonic function F as 



5(x) = F[<A(x)], 



(15) 



with i;^(x) being a Gaussian-distributed random field with vanishing mean and unit vari- 
ance, so that one-point PDF of 5 reads 



P[5]d5 
In this case we find 



1 



27r|F'[F-i(5)]| ''''^ I 2 



F-\d) 



d6. 



(16) 



r/i(x) = F'[</)(x)]</>,i(x), 

Ci,-(x) = F"[0(x)]</.,,(x)<A,,-(x)+F'[</.(x)]</.,,,(x). 



(17) 
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Therefore G(z/) is calculated as 



{[F"ml+F'm,^^] [F"(0)</.f2 + F'(</')</',22] - [F"(0)</,,i,/,,2 + F'(</,)0,i2]'} 



1 / 1 



2 \|F'(0)| 



<5d 



\F'm \F'{4>)\ 



(27r)2 \^ 



3/2 



exp 



1 



F-^{va) 



(18) 



where we have made use of the assumption that F{(f)) is monotonic with nonvanishing 

F'{^). 

Thus the shape of the genus curve is obtained from that of the Gaussian distribution 
simply replacing v by F~^{ya). This is as expected because mapping in terms of a 
monotonic function does not change the shape of the density profile. Using the relations 
cj2 = (F(0)2)^ and erf = {F' {(jyf') ^{{V (f))^) ^ for isotropic spaces, the overall factor can 
be rewritten to yield the final result: 



= (m. 



3/2 



\{F'{ci>y 

In the particular case of the lognormal distribution, the function F is defined by 

1 



(19) 



F[0(x)] = 



: exp 



^ln(l + a2)0(x) 



- 1. 



(20) 



We therefore find the PDF in the lognormal distribution as 

{ln[(l + ;/(T)v/TT^]} 



^'ln(i^) 



(1 + l^a)^/W^T+a^ 



5^ exp 



21n(l+a2) 



(21) 



and the genus curve in the lognormal distribution as 
1 
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(2^)^ [3(l + C72)ln(l + C72)]=^/2 



X exp 



V 



{in [(1 + /yrT)v/TT^]} \ ( {in [(1 + ;/(t) v/TT^] } 



21n(l + (72 



ln(l + a2) 



(22) 



One can easily check that the above expression reduces to the Gaussian formula in the 
limit a — > 0. 
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4 The Case Density Fields Depends on Multiple Gaussian 
Fields 

4.1 General consideration 

We now extend the above analysis to the case statistical properties of the density field 
is characterized by a number of independent Gaussian fields q;"(x) through a function / 
as 

(5(x)=/[ai(x),a2(x),...,a"(x)]. (23) 

Here we assume that q;"'s are mutually independent and all have vanishing mean and 
unit variance. Then introducing new variables Pf = a"j and = a"jj we find 



da"" 



(24) 



Here and below summation over repeated Latin indices a,b,c,...,h is implicitly assumed. 
Then the genus number is formally written as 



Giiy) = -^{SB{f-iycT)6Mf3t)W9Pi)\Ml\ 
Here the two-body correlations are given by 



(25) 



2 

^2 xab/ 



where a\ and (72 are defined by 

al^iiVa^f), and ^ ((VV)^). 



(26) 



(27) 



Introducing new variables a;^- = 7?- + 5^ q;"/3 we find a", (3\, and are random- 
Gaussian variables whose correlations are totally decoupled from each other: 



raft 



' - 2 -"4 \ "2 



(28) 



The above fact greatly simplifies the subsequent averaging procedures, because we can 
average over (jf^-, (3^, Pf, and P2 turn independently as follows. 
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First, in 



G{u) 



(29) 



c pd 



averaging over ujfj can be readily done using (28). The expression in the curly bracket 
should be replaced by 

{■■■}). = - Y/aA('^"'-a"a'')-^/,,(/??/3j+/?f/52^)/,a^+Afe/,,(/3?/3?/32^/32"^ 

(30) 

Next /Jf— average is performed noting that the linear combination fhPs is also a Gaussian 
with a vanishing mean and the variance 



{{Mlf)p, = "ff, with r ^ 



a 



to yield 



(31) 



(32) 



We assume only nonvanishing / contributes to the final result. The averaging over /3f 
or is more involved but it can be done utilizing the fact that Uj = fa(3f, /3f , and 
(p 7^ q) constitute a trivariate Gaussian distribution with the correlation matrix 



/ fp fq 

fp 1 
fq 1 



(33) 



After some straight forward calculations we find 



2^/ai' 



3 ^ ^, f ^Pi fpfq -2 



'1' 



(34) 
(35) 



2nfa, 3/2 

including the case with p = q- We thus obtain the expression for G(i^) leaving only 
a— average as 

fabfafb 



P 



fca' + {fa 



fabfab ~„ {fabfafb fee facfbefafb) 
f^ 



(36) 



This average cannot be calculated until we specify the function f{a^). We therefore 
move on to a specific example of the chi-square distribution in the next subsection. 
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4.2 Chi-square distribution 

In chi-square matter distribution the density field, p(x), is given by 

_ n 

pW = ^E[«'W]'' (37) 

p=i 

with J) being the mean density. The density contrast reads 



5{^) = P^^ = l±[a''{^)f-l. (38) 
Hence we find 

a' = {6{^f) = l. (39) 
n 

Although n is a positive integer by definition, we can perform analytic continuation to 
arbitrary positive number and replace n by using (pg|) in what follows. 
Thus in the chi-square distribution the function /(a) is given by 

/(a) = ya'^a'^ (40) 
with fa = cr^a" and fab = a'^S^^. We therefore find from (|36|), 

G{u) = --d^/6n{f{d)-i^a)^\a^a^ + a^a^-4a^a^ + 2a''-6a^ + 4]) 
(ovr)'^/^ \ (T^a L J / ^ 

' ^-.^-^^Vch(.)-Gch(.). (41) 



(37r)3/2 ^/rrpi;^ V 2 
Here -Pch('^) is the one-point PDF of 6{x)/a as a function of the threshold z^, which is 



calculated as 



Pcni'^) = {5Ty{f{a)/a-v)) 



d^a" _£^2/2 /a. 2 1 

-e ' Od —a V 



(27r)"/2 a 
(1 1^0-)'""""^ / l-Fl/cr 



a? = 3(r?2) = 3(/,/??) = a^a?^ = 2a^al (43) 



Using the relation 



we finally obtain 



(44) 
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Since ( p2[ ) reduces to the Gaussian distribution for a — > also coincides with (^) 
in this hmit. 

Note that (^) indicates that one-point PDF of the density field obeys the negative 
binomial distribution 

The negative binomial distribution is often used to fit the data of counts-in-cells analysis, 
namely, as a model of a PDF of smoothed density field. A number of observational 
analyses have shown that it fits the counts-in-cells histogram well for some volume-limited 
redshift samples (Gaztafiaga & Yokoyama 1993; Bouchet et al. 1993). We also note that 
this distribution has a hierarchical property of higher-order connected moments. That 
is, A^— th order cumulant kn is given by kn = (N — 1)1{k2)^^^ = {N — 1)!(t^^~^. 

On the other hand, in spite of the above similarity, the chi-square distribution and 
the negative binomial distribution are essentially different from each other once spatial 
dependence is taken into account. To illustrate it let us consider two- and three-point 
correlation functions, ^ and respectively, in the chi-square distribution. The former is 
calculated as 

e(x,y) ^ (5(x)5(y)) = ^(a«(x)a'^(x)a^(y)a^(y)) - 1 = a'w\x,y), (46) 

where (a"(x)a^(y)) = 6"'^w{'K,'y). Similarly, the latter is given by 

C(x,y,z) ^ (<5(x)5(y)<5(z)) 

= 2cr'^ii;(x,y)'u;(y,z)w(z,x) (47) 
= 2a[e(x,y)e(y,z)e(z,x)]i/2^ 

which is very different from a hierarchical form. On the other hand, if the negative 
binomial distribution would fit the counts-in-cells analysis for any shape and size of the 
sampling cell, we would expect the following relation to hold 

C(x, y, z) = ^ [C(y, z)C(z, x) + C(z, x)e(x, y) + ^(x, y)e(y, z)] , (48) 

with the observed slope (Peebles & Groth 1975; Groth & Peebles 1977) of the power-law 
two-point correlation function (Gaztahaga & Yokoyama 1993). Thus the two distribution 
should be distinguished from each other even if one-point PDF has exactly the same form. 

5 Comparison with Numerical Simulations 

We measure the genus of the four data sets from cosmological A^-body simulations with 
random-Gaussian initial conditions, kindly provided by T. Suginohara and Y. Suto. 
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Three models are evolved in the Einstcin-de Sitter universe with the scale-free initial 
fluctuation spectra (at expansion factor a = 1.0): 

P{k) a A;" (n = -1, 0, and 1). (49) 

The last model corresponds to a spatially-flat low-density cold dark matter (LCDM) 
model. In this specific example, we assume = 0.2, Aq = 0.8, and h = 1.0 (Suginohara 
&: Suto 1991). The amplitude of the power spectrum in the LCDM model at a = 6 
is normalized so that the top-hat smoothed rms mass fluctuation is unity at Sh^^Mpc. 
In fact this LCDM model can be regarded to represent a specific example of the most 
successful cosmological scenarios so far (e.g., Suto 1993). All models are evolved with 
a hierarchical tree code implementing the fully periodic boundary condition in a cubic 
volume of L^. The physical comoving size of the computational box in the LCDM model 
is L = 100/i~^Mpc. The number of particles employed in the simulations is N = 64^, 
and the gravitational softening length is Cg = L/1280 in comoving. Further details of the 
simulation models and other extensive analyses are described in Hernquist, Bouchet &; 
Suto (1991), Suginohara et al. (1991), Suginohara & Suto (1991), Suto (1993), Matsubara 
& Suto (1994), and Suto & Matsubara (1994). 

The computation of the genus from the particle data is performed using the code 
kindly provided by David Weinberg (Weinberg 1988; Gott et al. 1989). In short the 
procedure goes as follows; (i) the computational box is divided into N^{= 128^) cubes, 
and the density Pg{r) at the center of each cell is computed using Cloud-In-Cell density 
assignment, (ii) the Fourier-transform: 

is convolved with the Gaussian filter, and transformed back to define a smoothed density 
of each cell (with the filtering length Rf): 

Ps{^\ Rf) = ^J Poi^) exp ( ^ - ik • r j d'k. (51) 

(iii) the rms amplitude of the density fluctuations is computed directly from the smoothed 
density: 

aiRf) ^ ^{{ps/p-iy), (52) 

where p is the mean density of the particles, (iv) The isodensity surface of the critical 
density: 

Pc=[l + i^c7{Rf)] p (53) 
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is approximated by the boundary surface of the high-density {ps > Pc) and low-density 
{ps < Pc) cells, (v) Then the genus of the surface is computed by summing up the angle 
deficit D{i,j, k) at the vertex of cell (i, j, k): 

5^('^) = -4^ E (54) 

i,j,k=l 

The way to compute D{i,j,k) is detailed in Gott et al. (1986). The genus curve G{i^) 
is defined to be the number of genus per unit volume as a function of the threshold u. 
(vi) We repeated the above procedure 50 times using the bootstrap resampling method 
(Ling, Frenk & Barrow 1986) in order to estimate the statistical errors of G{v). 

It should be noted that earlier papers (e.g., Gott et al. 1989; Rhoads, Gott, & 
Postman 1994; Vogeley et al. 1994) defined the density threshold v of genus curves so 
that the volume fraction on the high-density region of the isodensity surface is equal to 

/= 1 r e-'^'^dt. (55) 



2tt 

Adopting this method, the lognormal model and random-Gaussian model would not be 
distinguished from each other (Coles & Jones 1991). Since we intend to distinguish them 
so as to see whether or not the lognormal model can fit the genus number, we adopt the 
straightforward definition 6 = va the density threshold throughout this paper. 

To obtain the normalized genus curve G(z^)/G(0), we follow the method developed 
in Matsubara & Suto (1995). In practice, we first compute G{v) at 51 bins (in equal 
interval) for —3 < < 3. Then we estimate the amplitude of G(0) by x^-fitting the 7 
data points around = to the lognormal formula (|2^ ) so that thus computed value of 
G(0) is less affected by the statistical fluctuation at one data point. The result is almost 
insensitive to which fitting formula we use in estimating G(0). 

The PDFs, P{v)^ which have been similarly computed with 50 bootstrap resampling 
errors, and normalized genus curves, G(i^)/G(0), are plotted in Figures 1 to 3 for power- 
law models with n = —1, 0, and 1, respectively. We select three different sets of the 
expansion factor a (= 1 at the initial epoch) and the filtering length Rf for each model so 
that the resulting u{Rf) covers from weakly to fully nonlinear regimes. The upper panels 
show the PDF and lower panels show the genus curves. The results of random Gaussian 
distribution, (|5|) and @j, are plotted in dotted curves, those of second-order perturbation 
theory, (|l|) and (|) with (|l|) in dashed curves, lognormal formulae, (|2l| ) and (P2|), in 
solid curves, chi-square formulae, (|4^ ) and (^), in dot-dashed curves. Symbols indicate 
the results of A^-body simulations. The curves for random Gaussian and second-order 
perturbation theory using Edgeworth series are not guaranteed that the genus for the 
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negative density vanish. We have forced these theoretical curves to be zero in the relevant 
regions in the plot. 

The comparison of genus curves between the A^-body results and the theories based 
on the Edgeworth series are described in Matsubara & Suto (1995) in detail, so we 
do not repeat the detailed argument here. In short, these two curves agree well for 
— 0.2< v(T< 0.4 where perturbation theory is expected to be valid, but the extrapolation 
of the second-order formula beyond this regime does not work. 

As for PDFs, the lognormal model fits the simulation results fairly well from weakly 
to fully nonlinear regimes. The degree of agreement, however, depends on the initial 
power spectrum. It works best for n = model. Note that we use the Gaussian window. 
Usually, the lognormal model is applied to the count-in-cells analysis which corresponds 
to the top-hat window. For the top-hat window, n = — 1 model fits the lognormal model 
better than n = model (Bernardeau &; Kofman 1995). Other non-Gaussian curves 
including chi-square distribution fit simulation data well on weakly nonlinear regimes, 
while they deviate considerably from simulation results on fully nonlinear regimes. 

For the genus, the lognormal model is also the best among the non-Gaussian dis- 
tributions considered here. The degree of agreement between the lognormal model and 
simulation results also depends on the initial power spectrum. The fitting is the best 
for n = — 1 model. For other models, the lognormal model does not fit well in positive- 
threshold regions. For all the simulation results, the agreement of the lognormal model 
happen to be better than weakly nonlinear formula for low-density regions. This may 
be partly because the weakly nonlinear formula cannot naturally take into account the 
positivity of the density, but the lognormal model achieves it by construction. On the 
other hand, the chi-square model does not agree well with any simulation results except 
for the initial linear regime, even if it assures the positivity of the density. 

In Figure 4 are plotted the PDFs and normalized genus curves for LCDM model for 
an example of a realistic cosmological scenario. The smoothing length R is 4/i~^Mpc. 
If galaxies trace mass, a = 6 corresponds to the present epoch [z = 0). Thus a = 4 
and 5 correspond to z = 0.5 and 0.2, respectively. We find our simulation results of 
LCDM model and the lognormal model for Gaussian smoothed PDFs agrees fairly well 
except for some differences of the peak height. See also Ueda & Yokoyama (1995) for 
top-hat smoothed PDF. For the normalized genus curves, the lognormal model fits the 
simulation results well. 
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6 Conclusions and Discussion 



In the present paper, we have derived theoretical prediction of the genus statistics for 
some non-Gaussian distributions. We have generally calculated the expectation value of 
genus number for non-Gaussian fields that are given by nontrivial functions of Gaussian 
random fields. Two specific fields of this category, the lognormal distribution and the 
chi-square distribution were investigated in detail. 

As is seen in the figures the one-point PDF of the smoothed density field is fitted by 
the lognormal model fairly well for all the simulated models adopted here, i.e., power- 
low models with n = —1,0,1 and the LCDM model. The lognormal formula for the 
genus curve also works well to fit all the simulation data in the low density regions. 
But considerable deviation is observed in the positive-threshold regions in the nonlinear 
regime of the power-law models with n = and 1. On the other hand, the formula fits 
the entire regions of n = —1 and the LCDM models fairly well. Thus the genus statistics 
arc more appropriate than one-point PDF to distinguish between various initial power 
spectra for Gaussian smoothed field. This is because the former depends on the spatial 
derivatives of density field, too. 

Several arguments exist to explain the validity of the lognormal model as a statistical 
distribution in nonlinear regimes. Coles & Jones (1991) argues that the lognormal distri- 
bution is obtained from the continuity equation in the nonlinear regime but with linear 
or Gaussian velocity fluctuations, so that it may be adopted as a model to describe statis- 
tics in the weakly nonlinear regime. On the other hand, Bernardeau & Kofman (1995) 
claims its successful fit to the PDF of CDM type simulations is just a coincidence due to 
the particular shape of the CDM power spectrum based on the top-hat smoothing. Since 
LCDM power spectrum has a similar structure to n = — 1 power law in the interested 
scales, our results of the genus curves are consistent with these arguments. 

Observationally, the negative binomial distribution, which has the same one-point 
PDF as the chi-square distribution, has been shown to fit the counts-in-cells of various 
redshift samples well (Gaztaiiaga & Yokoyama 1993; Bouchet et al. 1993). However, 
our results indicate that chi-square formulae reproduce neither one-point PDF nor genus 
curve of the simulations. Does that imply these A^-body simulations have nothing to do 
with the real universe? Not necessarily, because here we are extracting information on 
relatively small length scales R ^ 4/j.^^Mpc using all the 64^ particles in the simulations, 
while the observational counts analysis has been performed on larger scales using volume- 
limited samples with much smaller number density of galaxies. The effect of sparse 
sampling in the counts analysis has been analyzed by Ueda & Yokoyama (1995) for the 
LCDM model and it has been shown that the chi-square or the negative binomial PDF 
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does not fit the data well if we use all the particles in the simulation but it fits well if 
we use sparsely-sampled data with a similar number density to that of volume-limited 
samples currently available. 

Thus our results do not rule out the LCDM model. On the contrary, it remains one 
of the most promising models of our universe (e.g., Suto 1993). To test the LCDM model 
further, we can use the genus curve by examining if it is fitted by the lognormal formula. 
Although the presently available redshift data are not statistically significant enough 
(e.g., Vogeley et al. 1994) to extract a specific conclusion, we can reasonably expect the 
statistical significance of observed genus curve will improve rapidly in the near future. 

Now that we have obtained a number of theoretical formulae for the genus curve we 
can make use of it not only to test the Gaussianity of the primordial fluctuations but 
also to discriminate between various models of structure formation. 
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FIGURE CAPTIONS 



Figure 1 The PDFs (upper panels) and normalized genus curves (lower panels) from 
the A^-body simulation data for n = — 1 power-law model (J^o = Ij = 0) are 
plotted by open circles. Three different sets of the expansion factor a (= 1 at 
the initial epoch) are selected and the Gaussian window function with the filtering 
length Rf = L/25 is used. The values of the expansion factor and the resulting 
variance a are indicated in the figure. The theoretical prediction for the random 
Gaussian field is plotted by dotted lines, lognormal model by solid lines, chi-square 
model by dot-dashed lines, second-order perturbation theory by dashed lines. 

Figure 2 Same as Fig. |6| for n = power-law model. 

Figure 3 Same as Fig. |6| for n = 1 power-law model. 

Figure 4 Same as Fig. |6| for low-density cold dark matter model (JIq = 0.2, Aq = 
0.8, h = 1.0) The adopted Gaussian filtering length Rj corresponds to 4/i~^Mpc 
(comoving). The top-hat smoothed rms mass fluctuation at a = 6 is unity at 
8/i-^Mpc. 
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This figure "figl-l.png" is available in "png" format from: 
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